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The problem of competitive nucleation in the framework of Probabilistic Cellular Automata is 
studied from the dynamical point of view. The dependence of the metastability scenario on the self- 
interaction is discussed. An intermediate metastable phase, made of two flip-flopping chessboard 
configurations, shows up depending on the ratio between the magnetic field and the self-interaction. 
A behavior similar to the one of the stochastic Blume-Capel model with Glauber dynamics is found. 
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Metastable states are common in nature; they show 
up in connection with first order phase transitions. Well 
known examples are super-cooled and super-heated liq- 
uids. Their statistical mechanics description revealed to 
be a challenging task. An approach based on equilibrium 
states has been developed via analytic continuation tech- 
niques [l[ and via the introduction of equilibrium systems 
on suitably restricted sets of configurations [1, H, Q ■ The 
purely dynamical point of view, dating back to Ref. 
has been developed via the pathwise technique Q and 
the potential theoretical approach ■ 

We shall stick to the dynamical description to inves- 
tigate competing metastable states. This situation arise 
in many physical processes, such as the crystalization of 
proteins [j| Q and their approach to equilibrium [Toj | . 
The extreme situation is represented by the glasses, in 
which the presence of a huge number of minima of the 
energy landscape prevents the system from reaching the 
equilibrium The study of these systems is difficult, 
since the minima of the energy and the decay pathways 
between them change when the control parameters are 
varied. It is then of interest the study of models in which 
a complete control of the variations induced on the en- 
ergy landscape by changes in the parameters is possible. 

In this perspective, the analysis of the Blume-Capel 
model in Ref. [H, EH and that of the Potts model in Ref. 
[IH are of great interest. In the Blume-Capel model the 
sites of the lattice can be either empty or occupied by 
a 1/2-spin particle. The interaction favors the presence 
of neighboring aligned spins; the chemical potential A 
controls the tendency to have particles or lacunas on the 
lattice and the magnetic field h, depending on its sign, 
favors either the pluses or the minuses. Depending on 
the parameters, in the zero temperature limit the stable 
state is the one with all the spins up (u) or all the spins 
down (d) or no particle at all (0). Let h, A > 0, so 
that the unique stable state is u, and set a = h/X. For 
a < 1 the transition from the metastable state d to u is 
achieved via a sequence of increasing plus square droplets 
in the sea of minuses. For 1 < a < 2 and h small, the 



transition from d to u is realized via increasing squared 
frames in which the internal pluses are separated by the 
external minuses by a frame of zeros large one. For a > 2 
and h small, the system started at d visits the state 
before reaching u; the transition from d to is achieved 
via increasing zero square droplets in the sea of minuses, 
while the transition from to u is realized via increasing 
plus square droplets in the sea of zeros. 

We study, here, metastability for a Probabilistic Cel- 
lular Automaton [15( with self-interaction k, focusing on 
the dependence of the metastability scenario on such a 
parameter. The model interpolates those studied in Ref. 



[16| (k = 0) and [I7|,ll8| (k = 1). For k = each spin in- 



teracts only with its nearest neighbors; for k = 1 the self- 
interaction has the same strength as the nearest neighbor 
coupling. In absence of self-interaction an intermediate 
metastable state shows up; it is proven that the inter- 
mediate state is visited during the transition from the 
metastable to the stable state. The role played by the 
intermediate state changes as the self-interaction k is 
varied. Quite surprisingly, results similar to those found 
in Ref. |12j for the Blume-Capel model are obtained. 

Consider the two-dimensional torus A = {0, . . . , L — 
l} 2 , with L even, endowed with the Euclidean metric; 
x, y G A are nearest neighbors iff their mutual distance 
is equal to 1. Associate a variable o~(x) = ±1 with each 
site x G A and let S = { — 1,+1} A be the configuration 
space. Let j3 > and k, h G [0, 1]. Consider the Markov 
0, 1, ... , on S with transition matrix 



chain cr n , with n 



p(a, rj) 



Up- 

x6A 



(r)(x)) Vcr,?7e5 



(1) 



where, for x G A and a G S, p x .a(-) is the probabil- 
ity measure on { — 1,4-1} defined as p x .a{s) = 1/[1 + 
cxp{-2/3s(S' [T (x) +h)}} with s G {-1,4-1} and S a (x) = 
YlyeA K ( x - V) a (y) wner e K(x - y) is if \x - y\ > 2, 1 
if \x — y\ = 1, and k if \x — y\ = 0. The probability p XiCr (s) 
for the spin a(x) to be equal to s depends only on the 
values of the spins of a in the five site cross centered at x. 
The metastable behavior of model ([T]) has been studied 
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in Ref. [16( for k = and in Ref. [17|, for k = 1. 

The Markov chain {T]) is a probabilistic cellular au- 
tomata; the chain <t„, with n = 0,1,..., updates all 
the spins simultaneously and independently at any time. 
The chain is reversible, see Ref. (lBj . with respect to 
the Gibbs measure /i(a) = exp{— /3H(a)}/Z with Z = 
Tories ex P{~P H (v)} and 

H(a) = -hJ2 <t(x)~ Yl lo S cosh $ ( 5 -(- T ) + h )l ( 2 ) 



that is detailed balance p(a, rf) e~^- ff ( <T ) = p(rj, a) e~P H ^ 
holds for a,r] £ S; hence, /j, is stationary. We refer to 1//3 
as to the temperature and to /i as to the magnetic field; 
the interaction is short range and it is possible to extract 
the potentials as described in Ref. (l8j . 

Although the dynamics is reversible w.r.t. the Gibbs 
measure associated to the Hamiltonian <j2j) , the probabil- 
ity p(o~, 77) cannot be expressed in terms of H (a) — H(rj), 
as usually happens for Glauber dynamics. Given a, 77 G 
<S, we define the energy cost 



AM) 



lim 

/3->oo 



logp(cr, 7?) 



= J2 ^\Sa(x) + h\ (3) 



r,(x)[S o -(x) + h]<0 



Note that A(cr, 77) > and A(cr, 77) is not necessarily equal 
to A(t7, <t); it can be proven, see [l?], Section 2.6], that 

e -PM",v)-fhW < ^) < e -PM*,v)+fhW (4) 

with 7(/3) — > in the zero temperature limit (3 — » 00. 
Hence, A can be interpreted as the cost of the transition 
from (j to 77 and plays the role that, in the context of 
Glauber dynamics, is played by the difference of energy. 

To pose the problem of metastability it is necessary to 
understand the structure of the ground states; since the 
Hamiltonian depends on /?, their definition deserves some 
thinking. The ground states are those configurations on 
which the Gibbs measure /i concentrates when f3 — > 00; 
hence, they can be defined as the minima of the energy 

E{a) = lim H{a) = -/iV a(x) - V \S a (x) + h\ (5) 
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For X C S, wc set E(X) = mm aeX E(a). For h > the 
configuration u, with u(x) = +1 for x G A, is the unique 
ground state, indeed each site contributes to the energy 
with —h — (4 + k + h). For h = 0, the ground states are 
the configurations such that all the sites contribute to the 
sum ([5]) with 4 + k. Hence, for k £ (0, 1], the sole ground 
states are the configurations u and d, with d(x) = — 1 
for x G A. For n = 0, the configurations c e , c° G S 
such that c c (x) = and c°(x) = (-1)^+^+1 

for x = (xi,X2) G A are ground states, as well. Notice 
that c e and c° are chessboard-like states with the pluses 
on the even and odd sub-lattices, respectively; we set 



c = {c c ,c°}. Since the side length L of the torus A is 
even, then £(c°) = E(c°) = E(c). 

We study those energies as a function of k and h, re- 
calling that periodic boundary conditions are considered. 
Wehave£(u) = -L 2 (4+n+2h), E(d) = -L 2 (4+n-2h), 



and E(c) = -i 2 (4 



hence E(c) > E(d) > E(u) for 



< h< k < 1, E(c) = E(d) > E(u) for < h = k < 1, 
and E(d) > E(c) > E(u) for < K < h < 1. 

We can now pose the problem of metastability at fi- 
nite volume and temperature tending to zero (Friedlin- 
Wentzel regime). Following Ref. see also Ref. (l7l 
Appendix], given a sequence of configurations w = 
u>l, . . . , ui n , with n > 2, we define the energy height along 
the path to as $ w = max i=1) ... i | w |_ 1 [S(a;i) + A(a;i,Wj + i)]. 
Note that the definition does not depend on the direction 
in which the path u> is followed. More precisely, denoted 
by uj' the path uj n , w n _i, . . . , u>i, since 



E(a) + A(<7, 77) = £(77) + A(t ? , a) 



(6) 



for any cr, 77 G 5, it follows that $ w = $ W '; © is conse- 
quence of the detailed balance principle. Given A, A' C 
iS, we let the communication energy between A and A' 
be the minimal energy height $ w over the set of paths w 
starting in A and ending in A'. For any a G 5, we let 
Xo- C 5 be the set of configurations with energy strictly 
below E(a) and V a — ^{a,I a ) — E(a) be the stability 
level of <7, that is the energy barrier that, starting from 
<7, must be overcome to reach the set of configurations 
with energy smaller than E(a); we set V a = 00 if T a = 0. 
We denote by <S 8 the set of global minima of the energy 
([5]); namely, the collection of the ground states, and sup- 
pose that the communication energy T = max ff6 5^ B V a 
is strictly positive. Finally, we define the set of metastable 
states S m = {77 G S : V v = T}. The set S m de- 
serves its name, since it is proven the following (see, 
e.g., Ref. [13, Theorem A. 2]): pick a G 5 m , consider 
the chain o~ n started at o~q = a, then the first hitting 
time t$ s = inf{£ > : a t G 5 s } to the ground states is 
a random variable with mean exponentially large in j3, 
that is 



lim - logE CT [7>] =r 

0^oc [j 



(7) 



with E CT the average on the trajectories started at a. 

In this regime the description of metastability is re- 
duced to the computation of S s , T, and 5 m . We choose 
the parameters of the model (fl]) in such a way that 
< h < 1, h ^ K, and 2/h, 2/(h - k), 2/(h + k), and 
(2 + k — h)/h are not integer. The configuration u is then 
the unique ground state, i.e., S s = {u}. Two candidates 
for metastability are d and c; to find S m , one should 
compute r and prove that either Vd or V c is equal to Y. 
This is a difficult task, indeed all the paths u> connecting 
d and c to u must be taken into account and the related 
energy heights $ w computed. Since at each time step 
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all the spins of the lattice can be updated, the structure 
of the trajectories is highly complicated. This is why 
the study of the energy landscape of probabilistic cellu- 
lar automata is very difficult [ljj, Theorem 2.3]; such a 
task is simpler for serial Glauber dynamics, where a sort 
of general approach can be developed @, Section 7.6]. 

We develop an heuristic argument to compute T. Re- 
call ([3]) and note that k and h have been chosen so that 
S a (x) + h ^ 0. Thus, it follows that, given ue5, there 
exists a unique rj G S such that A(cr, n) = 0; the config- 
uration rj is such that rj(x){S a (x) + h] > for all x G A 
and is the unique configuration to which the system can 
jump, starting from cr, with probability tending to one 
in the limit j3 — > oo (sec ((!])). We say that a G S is a lo- 
cal minimum of the energy iff A(cr, a) = 0; starting from 
a local minimum, transitions to different configurations 
have strictly positive energy cost and thus happen with 
negligible probability in the zero temperature limit. It is 
immediate that d and u are local minima of the energy, 
while c G and c° are not, indeed c e (x) [S c ° (x) + h] < and 
c°(x)[S c o(x) + h] < for all x G A. We also have that 
A(c c ,c°) = A(c°,c c ) = 0, hence at very low tempera- 
ture, the system started in c° is trapped in a continuous 
flip-flop between c° and c e . A peculiarity of parallel dy- 
namics is the existence of pairs a, rj G S in which the 
chain is trapped since A(cr, n) = A(?], a) = 0; the proba- 
bility to exit such a pair is exponentially small in (3. 

We characterize, now, the local minima and the trap- 
ping pairs. For what concerns the local minima, we con- 
sider a configuration cr and study the sign of S (T (x) + h. 
Suppose, first, h < n and recall K < 1; the sign of 
S a (x) + h equals the sign of the majority of the spins 
in the five site cross centered at x. Hence, cr is a local 
minimum iff for each site x there exist at least two near- 
est neighbors such that the associated spins are equal to 
a(x). Suppose, now, h > k > 0; the sign of S a (x) + h 
is negative iff at least three among the spins associated 
to neighboring sites of x are minus. Hence, cr is a local 
minimum iff for each site x such that o~(x) = — 1 there 
exist at least three negative minus neighbors and for each 
site x such that <j(x) = +1 there exist at least two posi- 
tive neighbors. In conclusion, for h > k the local minima 
of the energy are those configurations in which all the 
pluses, if any, are precisely those associated with the sites 
inside a rectangle (plus-minus droplets). For h < k the 
local minima are all the configurations that can be drawn 
adding pluses to d so that each plus (rcsp. minus) has at 
least (resp. at most) two neighboring pluses. Plus-minus 
rectangular droplets arc local minima also in this case. 
For what concerns the trapping pairs, consider a configu- 
ration tr with a rectangle of chessboard plunged in the sea 
of minuses (chessboard-minus droplet) and let rj be the 
configuration obtained flipping all the spins associated 
with sites in the chessboard rectangle. The configuration 
cr, r\ form a trapping pair only for h > k. Indeed, it is 
immediate to show that all the spins of the chessboard 



tend to flip, some thinking is necessary only for the mi- 
nus corners. Let x be the corner site with cr(x) = — 1, 
since S a (x) + h = —n + h, we have that S a (x) + h > 
for h > k and S a (x) + h < for h < k. Thus, the spin 
tends to flip in the former case and not in the latter. 

The local minima and the trapping pairs can be used 
to construct the optimal paths connecting d and c to the 
ground state u. We distinguish two cases. 

Case h > k > 0. Although c c and c° arc not local 
minima of the energy, the system started in c is trapped 
in a continuous flip-flop between c° and c c . This trap- 
ping persists even if a rectangle of pluses is inserted in 
the chessboard background (plus-chessboard droplet); a 
path from c to u can be constructed with a sequence of 
such droplets. The difference of energy between two plus- 
chessboard droplets with side lengths respectively given 
by £,m > 2 and £, m + 1 is equal to 4 — 2(k + h)£. It then 
follows that the energy of a such a droplet is increased 
by adding an £-\ong slice iff t > [2/(k + ft)J + 1 = A" 
([x\ denotes the largest integer smaller than the real x). 
The length A" is called the critical length. It is reasonable 
that the energy barrier V c is given by the difference of en- 
ergy between the smallest supercritical plus-chessboard 
droplet, i.e., the plus-chessboard square droplet with 
side length A", and the configuration c; by using ([5]) 
we get that such a difference of energy is equal [19( to 

rs = 8/(« + /i). 

A path from d to u can be constructed with a sequence 
of plus-minus droplets. By using we get that the dif- 
ference of energy between two plus-minus droplets with 
side lengths respectively given by £, m > 2 and £, m + 1 
is 4(2 — hi). It then follows that the energy of a plus- 
minus droplet is increased by adding an £-long slice iff 
£> [2/h\ + 1 = Aj. The length A^ is the critical length 
for the plus-minus droplets; by using ((5|) we get that the 
difference of energy between the smallest supercritical 
plus minus droplet and d is equal to rjj = 16/ h. 

An alternative path from d to u can be constructed 
via a sequence of frames with the internal rectangle of 
pluses separated by the external minuses by a stripe of 
chessboard large one. These are peculiar trapping pairs 
in which the flip-flopping spins are those associated with 
the sites in the stripe of chessboard. We can prove that 
the difference of energy between two frames with internal 
(rectangle of pluses) side lengths respectively given by 
£, m > 2 and £, m + 1 is equal to 8 — A(h — k) — Ah£, 
so that the critical length for those frames is given by 
X f d = [(2 — h + K)/h\ + 1 and the difference of energy 
between the smallest supercritical frame and d is equal 
to I* = 16[1- (h-n)/2} 2 /h. 

A path from d to c can be constructed with a sequence 
of chessboard-minus droplets. By using ([5]) we get that 
the difference of energy between two chessboard-minus 
droplets with side lengths respectively given by £, m > 2 
and £,m+l is equal to 4 — 2(h — k)£. It then follows that 
the energy of a chessboard-minus droplet is increased by 
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FIG. 1: The time unit is the time step of the chain. Solid lines 
(from the left to the right) represent the magnetization of the 
runs («, 0) = (0.15, 0.55), (0.4, 0.5), (0, 025, 0.7). Dashed lines 
represent the absolute value of the staggered magnetization; 
the non-null curve is found for (k,/3) = (0.025,0.7). 

adding an £-long slice iff £ > [2/(h - k)J + 1 = \%. 
The length A d is the critical length for the chessboard- 
minus droplets; the energy difference of energy between 
the smallest supercritical chessboard-minus droplet and 
d is equal to r ^ = 8/(h — k). 

Note that Tjj < T d for h, K small. Moreover, let a = 
h/n and remark that, provided the magnetic field h is 
chosen small enough as a function of a. T d < T d for 
a > 2 and r 5 > T f d for 1 < a < 2. Hence, for a > 2 



we get Vd = T d , that is the chain escapes from d and 
reaches the state c in a time that can be estimated as in 
([7]) with r = rjj. Starting from c the chain will reach u 
by overcoming the energy barrier V c = T" < Va- Note 
that V c = Vd in the limiting case k = , hence both c 



and d are metastable states (results in [16j are recovered). 
For l<a<2,Vd=r^, that is the chain escapes from 
d and reaches the state u via a sequence of increasing 
frames in a time estimated as in (JT]) with r = T f d . 

Case h < k < 1. By paying the smallest energy cost 
any local minimum can be transformed in a configura- 
tion with the pluses forming well separated rectangles 
(see 18[); hence, the most relevant local minima are 



the plus rectangular droplets. As noted above, for this 
choice of the parameters the system cannot be trapped 
in chessboard-minus droplets. Thus, the energy barrier 
Vd is given by the energy T d of the smallest supercritical 
plus droplet. As before, we also have V c = T". Since 
V c < Vd, we have that d is the unique metastable state, 
the communication energy is T = T d , the tunneling time 
is exp{/?rjj} in the sense (0, and the zero temperature 
limit transition from the metastable state d to the sta- 
ble state u is achieved via the nucleation of a plus-minus 
square droplet with side length X d . For k = 1 the results 
proven in [171 ] are recovered. 

The mctastability scenario depends on the ratio be- 



tween the magnetic field and the self-interaction. For 
K = the two states d and c are both metastable. 
For a > 2 and h small, c is crucial, although not 
metastable, since it is visited during the transition from 
the metastable state d to the stable state u. For 2 > 
a > 1 and h small, the chessboard configuration plays 
no role at all and the exit from the metastable d state 
is achieved via the direct formation of the plus phase 
via a sequence of increasing frames. For 1 > a, the exit 
from the metastable d state is achieved via the direct 
formation of the plus phase via a sequence of increasing 
plus-minus droplets. The scenario is very similar to the 
one proven in Rcf. fl2| for the Blume-Capel model with 
Glauber (serial) dynamics; the role of the chemical po- 
tential A is played here by the self-interaction k. This be- 
havior has been tested at finite temperature via a Monte 
Carlo simulation [2(j. We have considered L = 1000, 
h = 0.2, and run the chain for («;,/?) = (0.025,0.7), 
(0.15,0.55), (0,4,0.5). By measuring the staggered and 
the usual magnetization, we point out that the system 
visits c before reaching u only in the run k = 0.025 and 
j3 = 0.7 (see Figure [I}, which is the only run with a > 2. 
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